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ABSTRACT 


This  thesis  consists  of  two  main  parts.  The  first  is  an 
organized  presentation  of  some  results  in  queuing  theory  that  are 
applicable  in  the  analysis  of  telephone  trunking  problems.  Here 
four  situations,  all  assuming  Poisson  input  and  exponential  holding 
times,  are  discussed.  These  are  the  case  of  the  f irst-come-first- 
served  single  channel  queue,  the  case  of  an  infinite  number  of 
servers,  the  case  of  a  finite  number  of  servers  where  one  waiting 
line  is  formed  on  a  f irst-come-f irst-served  basis  and  finally  the 
case  of  a  finite  number  of  servers  where  an  arriving  item  is  lost 
if  all  servers  are  busy. 

The  second  part  of  the  thesis  is  an  attempt  to  analyse 
some  actual  data  obtained  from  Alberta  Government  Telephones  on 
the  calls  terminating  at  the  699  77 —  Sherwood  Park  telephones 
from  9  A.M.  to  5  P.M.  over  about  an  eight  day  period.  Here  the 
interarrival  times  and  busy  times  are  calculated  and  the  results 
compared  with  the  exponential  assumptions  made  previously. 
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CHAPTER  I 


INTRODUCTION 


1.1  A  Discussion  of  Queuing  Theory 

When  people  arrive  at  a  facility  that  provides  some  necessary 
or  desirable  service,  whether  it  be  a  check-out  counter  in  a  grocery 
store  or  a  subway  station,  they  frequently  encounter  a  queue  or  waiting 
line.  The  study  of  these  waiting  lines,  queuing  theory,  is  a  mathemat¬ 
ical  descriptive  theory  that  attempts  to  formulate,  interpret  and 
predict  queuing  behavior  for  the  purpose  of  better  understanding  this 
bothersome,  but  seemingly  inevitable,  aspect  of  our  daily  lives  and  for 
the  sake  of  introducing  remedies . 

The  importance  of  studying  and  developing  such  a  theory  is 
clearly  evident  in  the  amount  of  time  an  individual  is  made  to  wait 
in  his  daily  activities.  If  one  is  not  properly  occupied  many 
queuing  situations  amount  to  an  appreciable  depreciation  of  the 
enjoyment  of  life.  But  by  anticipation  and  good  planning  the  agony 
and  waste  of  queues  could  be  minimized. 

Queues,  however,  are  an  unavoidable  part  of  our  modern  society. 
With  an  ever  increasing  population,  the  number  of  people  demanding  the 
same  services  will  continue  to  grow  and  unless  these  services  are 
increased,  which  in  some  cases  can  only  be  done  at  enormous  costs,  the 
amount  of  time  an  individual  spends  in  a  queue  in  his  lifetime  will 
increase  considerably.  Thus  it  becomes  evermore  necessary  to  understand 
and  anticipate  queues  and  their  subtleties. 
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However,  in  attempting  to  do  this  it  must  be  kept  in  mind  that 
most  queuing  systems  involve  human  beings.  Studying  causes  and  remedies 
of  queuing  problems  cannot  be  completely  divorced  from  consideration  of 
human  factors  and  their  influence  on  the  problem.  Nor  can  remedies  be 
applied  without  regard  to  the  fact  that  it  is  the  people  using  the  facil¬ 
ities  who  matter  in  the  final  analysis. 

1.2  The  Application  of  Queuing  Theory  to  Telephony 

It  is  interesting  and  perhaps  of  particular  importance  that 
the  origin  of  queuing  theory  is  to  be  found  in  telephone-network 
congestion  problems.  Some  basic  ideas  of  telephony  that  are  of 
particular  interest  to  the  remainder  of  this  presentation  are  described 
as  follows. 

Let  us  assume  that  we  have  a  set  of  devices  (telephones, 
trunks,  lines,  etc.)  that  are  all  accessible  to  a  group  of  sources 
(telephone  subscribers) .  When  a  caller  from  a  certain  group  wishes  to 
make  a  connection,  he  may  do  so  through  any  one  of  the  free  lines.  If 
all  lines  are  engaged  then  the  caller  either  waits  or  his  call  is  lost; 
that  is  he  receives  a  busy  signal.  In  either  case,  the  probability  of  a 
call  waiting  a  certain  time  or  getting  lost  is  of  interest. 

Much  of  the  early  work  on  telephone  theory  was  done  by 
A.K.  Erlang.1  If  there  is  a  limited  source  of  N  subscribers  of  whom 

^Brockmeyer,  E.,  H.L.  Halstrom  and  A.  Jensen:  "The  Life  and 
Works  of  A.K.  Erlang"  (all  of  Erlang's  papers  translated  into  English), 
Copenhagen  Telephone  Company,  Copenhagen,  1948. 
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n  are  occupied  making  calls,  Erlang  assumed  that  the  probability  of  a 
new  call  being  initiated  in  the  time  interval  (t,t+dt)  is  (N-n)Xdt 
where  X  is  the  call  rate  per  unit  time  of  a  caller.  If  no  considera¬ 
tion  is  given  to  the  number  of  occupied  independent  callers,  then  the 
probability  of  a  new  independent  call  being  initiated  is  Xdt,  leading 
to  the  exponential  distribution  Xe  ^  for  the  distribution  of 
intervals  between  successive  calls.  This  obviously  gives  a  Poisson 
distribution  (Xt)ne  ^t/n!  for  the  probability  of  exactly  n  calls 
arriving  during  time  t  with  a  mean  of  Xt  calls.  Similarly,  if 
1/y  is  the  average  duration  of  a  call,  the  probability  of  a  call 
terminating  during  (t,t+dt)  is  ydt  which  is  independent  of  the 
length  of  the  call  and  of  other  calls.  This  in  turn  implies  that  call 
durations  are  distributed  according  to  ye  ^Jt. 

Erlang  developed  telephone  traffic  theory,  using  these 
assumptions,  and.  obtained  results  for  the  probability  of  different 
numbers  of  calls  waiting  and  for  the  waiting  time  when  the  system  is  in 
equilibrium  and  the  probability  of  loss  for  the  loss  system.  He 
assumed  Poisson  inputs  from  unlimited  sources  and  either  exponential  or 
constant  holding  times. 

1.3  The  Relationship  to  Markov  Processes 

The  main  assumption  in  the  above  discussion  is  that  the 
interarrival  times  and  service  times  are  exponential.  This  simplifies 
the  analysis  of  the  situation  as  then  the  process  (X(t)},  where  X(t) 
is  the  number  of  calls  in  the  particular  system  we  are  studying  at 


time  t,  is  certainly  Markovian. 
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Even  if  only  one  of  the  interarrival  times  or  the  holding  times 
is  exponential  we  may  still  use  the  Markov  approach. 

If  the  interarrival  times  are  exponential,  the  (x(t)} 
process  is  considered  just  when  a  customer  finishes  service.  In  this 
way  we  generate  a  sequence  of  integer  values 

X(T1) ,X(T2) ,X(T  ) , . . . 


where  Tg ,T2 ,T3 ' * ' '  are  the  successive  times  of  completion  of  service. 
The  sequence  (X(T^)}  forms  a  discrete  time  process  that  (because  of 
the  Poisson  nature  of  the  input)  describes  a  Markov  chain. 

Similarly,  the  exponential  nature  of  the  service  times  will 
induce  a  Markov  chain  at  the  times  of  arrival  of  new  customers. 

However,  if  neither  the  interarrival  times  not  the  service 
times  are  exponential,  we  do  not  have  an  embedded  Markov  chain  and  the 
situation  is  much  more  difficult.  Very  little  has  been  written  about 
this  case,  except  perhaps  for  one  server  systems.  Some  examples  of 
contributions  to  this  problem  can  be  found  in  D.V.  Lindley's  paper, 

"The  Theory  of  Queues  with  a  Single  Server",  and  in  papers  by  S.O.  Rice 
and  V.E.  Benes  that  appear  in  the  Bell  System  Technical  Journals. 

1.4  Outline  of  the  Thesis 

The  ensuing  material  will  be  divided  into  two  parts.  In  the 
first  part  there  will  be  a  development  of  some  of  the  mathematical 
ideas  that  have  been  established  in  telephone  theory,  beginning  with  a 
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brief  discussion  of  Erlang's  original  model. 

Here  he  describes  the  situation  of  a  f irst-come-f irst-served 
single  channel  queue  with  a  Poisson  input  with  parameter  X  and 
exponential  holding  times  with  parameter  y  .  By  solving  a  set  of 
differential  equations  that  give  the  change  in  the  probability  of  n 
calls  in  the  system  at  any  time  t  with  respect  to  time  we  will  obtain 
an  expression  for  this  probability.  By  setting  the  time  derivates  equal 
to  zero  and  eliminating  time  from  these  differential  equations  it  will 
be  possible  to  describe  the  situation  where  the  system  eventually 
reaches  an  equilibrium  state  independent  of  time,  giving  the  probability 
that  n  calls  are  in  the  system  in  this  steady  state.  Then  we  will 
approach  the  question  of  waiting  times  in  order  to  obtain  the  probability 
of  an  arriving  call  being  delayed  longer  than  time  t.  In  this  way  it 
could  perhaps  be  determined  whether  it  would  be  profitable  for  a  customer 
to  wait  or  not. 

Following  this  discussion,  the  case  of  an  infinite  number  of 
channels  that  are  all  accessible  to  a  certain  group  of  subscribers,  and 
thus  a  situation  which  involves  no  loss  or  waiting,  will  be  considered. 
Again  the  probability  of  n  calls  in  the  system  at  any  time  t  and  in 
the  equilibrium  state  will  be  determined. 

Next  we  will  consider  the  situation  of  a  finite  number  of 
channels  that  are  all  accessible  to  a  certain  group  of  subscribers  where 
callers  will  form  a  single  f irst-come-f irst-served  waiting  line  if  all 
channels  are  occupied.  In  this  case,  it  is  quite  difficult  to  determine 
the  probability  of  n  items  in  the  system  at  time  t  from  the 
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differential  equations.  However,  once  this  has  been  accomplished,  a 
straightforward  discussion,  much  similar  to  that  of  Erlang's  model,  will 
be  given  to  determine  the  equilibrium  probabilities  and  the  waiting  times 
both  in  the  queue  and  in  the  system. 

Finally  the  situation  of  a  finite  number  of  channels  that  are 
all  accessible  to  a  certain  group  of  subscribers  where  each  caller  will 
hang  up  if  all  channels  are  occupied  will  be  considered.  Here  the 
differential  equations  describing  the  change  in  the  probability  of  n 
calls  in  the  system  at  any  time  t  with  respect  to  time  will  not  be 
solved  as  they  involve  only  a  slight  modification  in  those  of  the 
previous  case.  The  steady-state  probabilities  will  be  given  however, 
and  also  the  probability  of  a  call  being  lost  to  the  system  (Erlang's 
loss  formula)  will  be  obtained. 

In  the  second  part  of  this  presentation  an  attempt  will  be 
made  to  analyse  some  actual  data  obtained  from  Alberta  Government 
Telephones  on  the  calls  terminating  at  the  699  77 —  Sherwood  Park 
telephones  from  9  A.M.  to  5  P.M.  over  about  an  eight  day  period. 

Here  we  will  begin  with  an  explanation  of  the  data,  which  is  in  the 
form  of  a  pen  recorder  chart,  and  then  outline  the  steps  that  were 
followed  in  interpreting  this  data,  in  calculating  the  interarrival 
times  and  busy  times  and  in  estimating  the  necessary  parameters . 

Finally,  the  results  obtained  will  be  discussed  and  conclusions  drawn 
as  to  how  well  these  results  compare  with  the  exponential  assumptions 


made  previously. 
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CHAPTER  II 


SOME  RESULTS  IN  QUEUING  THEORY 

2.1  A  Discussion  of  Erlang's  Original  Model 

Turning  now  to  the  discussion  of  Erlang's  original  model,  let 
P  (t)  be  the  probability  of  n  calls  in  the  system  at  time  t*, 
that  is  one  customer  being  served  and  (n-1)  waiting.  We  immediately 
have  the  equations 

P  (t+At)  =  P  (t) [(1-AAt) (1-yAt)] 
n  n 

+  P  . (t) AAt  (1-yAt) 
n-1 

+  P  (t) yAt  (1-AAt)  +  o (At)  for  n  >  1  ,  (1) 

n+1  — 

P  (t+At)  =P  (t) (1-AAt)  +  P. (t) yAt ( 1-AAt)  +  o (At) 
o  o  1 

We  arrive  at  these  equations  by  observing  that  the  probability 
of  n  calls  in  the  system  at  time  t+At  is  equal  to  the  probability 
of  n  calls  at  time  t  and  no  arrivals  or  departures  during  At  or 
the  probability  of  (n-1)  calls  at  time  t  and  one  arrival  but  no 
departures  during  At  or  the  probability  of  (n  +  1)  calls  at  time  t 
and  one  departure  but  no  arrivals  during  At  or  the  probability  of  more 
than  one  change  during  At,  which  is  o(At). 

By  transposing  and  passing  to  the  limit  with  respect  to  At 


we  obtain  the  system: 
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dP  (t) 
n 

dt 


(X+y)  P  (t)  +  XP  (t)  +  yP  .  (t) 
n  n-1  n+1 


for  n  >  1 


(2) 


dP  (t) 
o 

dt 


=  -  XP  (t)  +  yP  (t) 
o  1 


r 


which  can  be  solved  to  obtain  an  expression  for  P  (t) . 

n 

If  we  are  able  to  determine  the  generating  function 


P(z,t)  =  \  Pn(t)zn,  |z|  <_  1  ,  (3) 

n=o 


the  probabilities  P  (t)  which  we  seek  are  the  coefficients  of  zR  for 

n 

every  n  .  However,  instead  of  determining  P(z,t)  we  will  obtain  an 

* 

expression  for  its  Laplace  transform  P  (z,s)  and  calculate  its 

coefficients,  observing  that  P  (s)  is  the  transform  of  P  (t) . 

n  n 

We  begin  by  multiplying  each  equation  of  (2)  by  zn+^~  and 
summing  over  n  ,  obtaining 


z 


3P (z , t) 
3 1 


(1-z) [ (y-Xz) P (z,t) ~yP  (t) ] . 

o 


(4) 


If  we  assume  P. (0)  =  1,  then  P(z,0)  =  z1. 

l 

Taking  the  Laplace  transform  of  this  expression  and  solving 
for  P*(z,s)  we  obtain 


y 

. 

1 

I 
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P* (z,s) 


z1+1-]i  (1-z)  P*  (s) 
_ o 

sz-  (1-z)  (y-Az) 


(5) 


Since  P  (z,s)  converges  inside  and  on  the  unit  circle  for 
Re(s)  >  0,  the  zeros  of  the  numerator  must  correspond  with  those  of  the 
denominator  in  this  domain. 

The  zeros  of  the  denominator  are 


a,  (s)  = 
k 


A+y+s± [ (A+y+s)  -4Ay] 
2A 


k  =  1,2 


(6) 


where  a^(s)  has  the  positive  sign  before  the  radical. 

In  order  to  determine  whether  these  zeros  are  inside  the  unit 
circle  we  will  apply  Rouche ' s  Theorem  to  the  denominator 


sz  -  (1-z)  (y-Az) 

letting  f(z) 

=  (A+y+s) z 

and 

g(z)  =  Az  +  y  . 

Clearly  for  | z | 

=  1,  |f (z)  |  > 

| g  (z)  |  and 

thus 

sz  -  (1-z)  (y-Az) 

has  the  same  number  of  zeros  inside  the  unit  circle  as  f (z) ,  namely 
one.  This  zero  must  be  a  (s)  since  |a  (s) |  <  |a  (s) |  and  there  is 
no  zero  on  | z |  =  1 . 

1  * 

Thus  the  numerator  z  -  y(l-z)P  (s)  must  vanish  for 

o 

z  =  ,  that  is 


oi't  -  y  (1-a  )P*  (s)  =  0, 
2  2  o 


implying 


. 

■ 


? 


. 
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i+1 


P*(s) 

o 


y (1-a2) 


(7) 


Substituting  this  into  the  expression  for  P  (z,s)  gives 


z1+1-[(l-z)aj+1/(l-a  )] 

P  (z,s)  =  - — - — - - - 

-A (z-a  ) (z-a  ) 


(8) 


as  a  +  a  =  — f^+s  f  a  a  =  y/A  and  s  =  -A (1-a.)  (1-a.)  . 
1  2  A  1  2  12 


When  the  numerator  is  multiplied  through  by  1  -  a  ,  it 


simplifies  and  factors  giving 


P  (z,s)  = 


,  ,  .  i  i-1  i.  .  .  .  i-1  i-2  i-1. 

(z-a  ) (z  +a  z  +. . .+a  ) -za  (z-a  ) (z  +a  z  +...+a  ) 

zli  Zl.  Z>  b 

Aa  (z-a  ) (1-z/a  ) (1-a  ) 


.  i  i-1  i,  ,  i-1  i-2  i-lx 

(z  +a2Z  + .  .  . +a2) -za2  (z  +ct2Z  t.-.+o^  ) 

Xa  (1-z/a  ) (1-a  ) 

.  .  i  i-1  i,  i+1 

(1-a  ) (z  ta^z  +...+a2)+a2 

Aa  (1-z/a  ) (l-a2) 


1  ,  i  i-1  is  r  /  /  % 

tst  (z  +v  +---+a2)  i  (z/ai> 

1  k=o 


a. 


i+1 


+ 


T - — - -  \  (z/a  ) 

Aa. (1-a.)  .  L  1 


1  2  k=o 


(9) 


00 

where  \  ( z/a  x )  k  = 


k=o 


as 


| z/a  |  <  1  - 


' 

)  ( J 


■ 


11 


is 


The  coefficient  of 


n 


z 


in  the  first  term  of  the  above  statement 


+ 


y/A  ,  (y/A) 


+ 


A  a 


n-i+1  ,  n-i+3  ,  n-i+5 


+  •  .  •  + 


(y/A) 


A  a 


A  a 


1 


A  a 


n+i+1  ' 
1 


n  >  1  r 


(10) 


...  i-n  ,  .  i-  (n-1) 

(y/A)  +  (y/A) _ + 

.  i-n+1  ..  i-n+3 

Aa  Aa. 


+  ■ 


(y/A) 


A  a 


i+n+1 


n  <  i , 


(11) 


and  the  coefficient  of  zU  in  the  second  term  of  the  above  statement  is 


a. 


i+1 


a 


Ac^+1(l-a2) 


A  a 


i+1 

2  2 
— —r  (1+a  +a  +. 
n+1  2  2 


.)  ,  as 


a2l  <  i  ' 


1  n+1 

=  j  (A/y) 


l  (vA)k  \  ■ 


k=n+i+2 


a 


1 


(12) 


Thus  for  n  >  i, 


P*(s)  =  1/A  [ - c— r  +  y/^ 

n  n-i+1  n-i+3 


-j  CO 

+  1+AL.+  (X/y)n+1 


a. 


a. 


a 


n+i+1 

1 


l  +/«k  -Tl 


k=n+i+2 


a. 


(13) 


and  thus  we  will  have  P  (t)  by  taking  the  inverse  transform  of  P  (s) , 

n  n 


that  is 


H  - 


■i  OS' 


,j 


. 


12 


P  (t) 
n 


1 

2Tri 


c+i°° 

/ 


c-i°° 


e  StP* (s)ds 
n 


(14) 


To  determine  this  we  observe  that 
Laplace  transform  of2 


s+/s2-4Ay  -v 
1  2  A  ] 


is  the 


(2A)Vv(2/aTo  Vt  1Iv(2/A^  t) 

=  v  ( /A/y  )Vt  ^  (2/Ay  t)  (15) 

where  I  (z)  is  the  modified  Bessel  function  of  the  first  kind  given 

v 

by 


I  (z)  =  l 
V  .  u 


(z/2) 


v+2k 


k=o 


k ! T ( v+k+1)  * 


(16) 


2 

,  .  ( A+y+s )  +  [  ( A+ii+s )  — 4Ay]  2  * 

Thus  as  a  (s)  =  - gr -  and  as  f  (s)  is 

1  2  A 

"fc 

the  Laplace  transform  of  f(t)  implies  that  f  (s+a)  is  the  Laplace 
transform  of  e  atf(t)  we  have 


-  ( A+y) t 

P  (t)  =  - - - -  [(A/IT  ) 

n  A 


n  1+1  (n-i+1)  t  1I  .  n  (2/Ay  t) 

n-i+1 


y/A  (A7y  ) 


n  1+3  (n-i+3)  t  XI  .  ,(2/a^  t)] 

n-i+3 


2  "Bateman  Manuscript  Project",  California  Institute  of 
Technology,  McGraw  Hill  Book  Co.,  Inc.,  New  York,  1954,  Vol.  I,  4.16,  (4). 


? , 


X.I- 
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+  ...  +  (y/A) 1  ( /A/y"  )n+1+1(n+i+l)t  XI  .  n  (2 Ay  t) 

n+i+1 


+  ( A/y ) n+1  l  (/jVA  )kkt  1I  (2/xTT  t)]  . 

k=n+i+2 


(17) 


2  . 

Substituting  the  well-known  relation  —  I  (z)  =  I  ,  (z) 

z  v  v  —  1 


I  .  (z)  one  has 
v+1 


~(A+y)t  _  ,  _  _ 

P  (t)  =- - - -  {(/A/y  )n  1  /Ay  [I  .  (2/x^t)-I  .  _(2/Ay  t) 

n  A  n-i  n-i+2 


+  (/x7m~  )n  1+1  Ay"  [I  .^(2/^7  t)-I  .  .(2/mT  t)] 

n-i+2  n-i+4 


+ 


+  (/aa7  )n  1+1  /xiT  [I  ( 2/XyT  t)-I  ( 2/Ay"  t)] 

n+i  n+i+2 


+  (A/y)  J  (/y/A  )k  /ATT  [Ik_1  (2/XiT  t)-I  (2/mT  t)  ]  } 

k=n+i+2 


But 


(18) 


(A/y) n+1  I  (/y?A  )k[ik_1(2/M7  t) -ik+1  (2/aTT  t)  ] 
k=n+i+2 


=  ( A/y )  n+1  [  ( /y/A  )  n+1+2ln+i+i  (2 /Ay  t)+iy7A 


(47a  )kIv(2/A^  t) 


k=n+i+2 


+  (/77a  )n+1+1I  .  _(2/Ay  t)  -  /A?y  I  (47*  )kI  (2/Ay  t) 

n+1+2  k=n+i+2  k 


‘ 
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(/Vu  )n  t)  +  (/A 7u  )n  1+1I  ^.^^(2/Xu  t) 

n+i+i  n+i+2 


+  (l-x/u)  U/p)  /Vp  y  (/ivx  )  I,  (2/Vp  t) 

k=n+i+2  k 


(19) 


Thus ,  for  n  >  i  , 


Pn(t)  =  e  U+y)t[(47A  )n  1In_i(2/A7  t) 


+  (47a  )1  n+1I  .  _  (2/Mi  tWl-A/u)  ( A/y) n  l  (>47a  )KI,  (2  47  t)  ]  . 

n+i+i  ,  v  „  K 

k=n+i+2 


When  n  <  i  we  have 


(20) 


P*(s) 

n 


1/A  [(uA).  +  - 


a. 


i-n+1 


a. 


i-n+3 


(y/A)  , ,  ,  x  n+1 

+  . .  .  +  .  „■  +  (A/y) 


a. 


i+n+1 


y  (p/A)k  i-  ] 


k=n+i+2 


a. 


Thus 


—  ( A+y )  t  ,  , 

P  (t)  =  - - - -  [  (y/A) 1  n  (  /aTp  ) 1  n  (i-n+1)  t  I.  _  (2 /Ay  t) 

n  A  i-n+1 


+  ...  +  (y/A  )1(/Vy  )X  n  (i+n+1)  t  Ii+n+1( 2y/^  t) 


+  (A/y)  +1  £  ( 4/A  )kkt  1Ik(2/Ay'  t)  ] 

k=n+i+2 


;  ) 


‘ 


15 


- ( A+y) t 

- 7 -  {(/aTu  )n  1  /Ay  [I.  (2/Ay  t)-I .  <2/Ayt)] 

A  l-n  i-n+2 

+  .  .  -  +  (/A/y  )  n  i+1  /Ty  [I  (2/xU  t)  -I  .  (2/Ay"  t)  ] 

n+i  n+i+2 

oo 

+  (A/y)n  1  £  (/y/A  )k  /Ay"  [I  (2/A*y  t)-I  (2/Ay  t)  ] 

k=n+i+2  k  1  k+i 

(21) 


which  agrees  with  the  result  for  n  >  i  as  I  =  I 

—  -n  n 

For  many  problems  time  dependence  is  not  the  typical  situation 

to  be  studied.  After  a  long  period  of  operation  the  system  acquires 

a  pattern  of  behaviour  which,  although  describable  in  terms  of 

probabilities,  does  not  depend  on  time.  Thus  as  t  00  this  situation 

may  be  expected  to  exist  for  many  systems.  The  question  of  whether 

such  p  =  lim  P  (t)  exist  is  a  main  concern  of  ergodic  analysis  where 
n  n 

t-*®° 

one  actually  looks  for  the  probabilities  p^  which  describe  the  steady 
or  equilibrium  state.  But  one  also  notes  that  the  requirement  that  the 
probabilities  P  (t)  do  not  change  with  time  is  the  definition  of  the 
steady  state. 

The  changing  of  P  (t)  with  respect  to  t  is  described  by  the 

derivative  P' (t) ,  and  the  steady  state  condition  requires  that 
n 

P'(t)  =  0.  Thus  essentially  we  have  two  ways  of  obtaining  the  steady- 
n 

state  probabilities: 


1.  from  P' (t)  =0  and 
n 


2. 


lim  P  (t)  . 
n 


' 


1 


X  < 


ic.r. 


16 


By  setting  the  time  derivatives  equal  to  zero  and  eliminating 
time  from  the  system  of  equations  (2)  we  obtain  after  transposing 


(X+p)pn  =  Apn_1  +  upn+1  ,  n  >  1  , 


Apo  =  UPl 


(22) 


From  here  we  obtain 


p  =  pRp  ,  where  p  =  A/y 
n  o 


(23) 


Since  T  p  =  1, 
L  n 


n=o 


OO 


n=o 


OO  CO  -Q 

r  n  r  n  ^ o 

=>  L  PQP  =  P0  1  P  =  JIZ  =  1  if  P  =  A/y  <  1 
n=o  n=o 


=>  P  =  pn(l-p)  if  p  =  A/y  <  1  . 
n 


(24) 


Thus  in  the  case  where  the  traffic  intensity  factor 
p  =  A/y  <  1  we  have  the  steady  state  probabilities  satisfying  a 
geometric  distribution.  In  the  case  where  A/y  >_  1  the  number  waiting 
in  the  system  would  become  infinite  with  time  and  no  steady  state  occurs. 

Continuing  in  the  case  A/y  <  1,  the  expected  number  in  the 


system  is  given  by 


<= 

' 
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00  00 

L  =  l  np^  =  (1-p)  l  npn 
n=o  n=o 


=  (1-P)P 


d_ 

dp 


I 

n=o 


n 

P 


P 

1-P  * 


(25) 


Fluctuations  in  the  number  waiting  can  in  fact  occur  in  the 
steady-state  and  this  can  best  be  seen  by  calculating  the  variance 


P  \  2 


l  n-L  p  =  £  n  p  -  (~~  ) 

L  n  L  n  1-p 

n=o  n=o 


OO 

/  -i  %  v  2  n  p  2 

=  (1-p)  l  n  p  -  (— ) 


n=o 


P  \  2 


(1-p>e^p^  l  p  -  <igr> 

n=o 


e_  +  2p  .  (_E_) 2 


1  p  (l-p)2  1  p 


r-+  <T^)2 

1-p  1-p 


=  L  +  L  • 


(26) 


The  expected  number  in  the  line  is  given  by 


L 

q 


V  (n-l)p  =  L 
L  „  n 


P 


1-P 


P 


2 


_P _ 

1-p  ' 


(27) 


? 


■ 
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The  logical  question  to  ask  at  this  point  is  how  long  will  the 
individuals  in  the  queue  have  to  wait?  To  answer  this  question  suppose 
that  an  arriving  customer  finds  n  customers  ahead  of  him.  His 
waiting  time  will  clearly  be  a  random  variable  that  is  the  sum  of  the 
service-time  random  variables  of  the  n  individuals  ahead.  We  know  that 
all  the  service-time  distributions  are  identical  exponential  distributions 
with  parameter  y  and  that  because  of  the  "forgetfullness"  property  of 
the  exponential  distribution,  the  remaining  service  time  of  the  customer 
in  service  satisfies  the  same  law  as  the  service  time  of  those  who  have 
not  yet  been  served. 

If  we  consider  T.  as  the  service  time  of  the  ith  customer 

1 

then 


P{  (n+l)st  customer  waits  for  a  time  <_  t} 

=  W  (t)  =  P{T  +T  +...+T  <  t}  . 

n  12  n  — 

The  characteristic  function  of  the  exponential  distribution  is 

( — - — )  and  thus  the  characteristic  function  of  the  distribution  of 
y-i0 

the  waiting  time  of  the  (n+l)st  customer  is  (  )  ,  thus  giving  a 

density  of 


,  , A  n-1  -yt 
y(yt)  e 

(n-1)  ! 


which  implies 


■  . 

I 
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W  (t) 
n 


t  .n-1  -yx 

/  MiMl  ,  ^  dx 


o 


(n-1)  ! 


(28) 


The  density  function  of  the  distribution  of  the  waiting  time 
T  of  any  customer  in  the  line  is  given  by 


w(t) 


I  Pnd-P) 

n=l 


.  n-1  -yt 

U(Ut)  e 

(n-1)  ! 


=  yp(l-p)e 


- (1-p) yt 


(29) 


and  thus  we  can  calculate  the  average  waiting  time  as 


/  yp  (1-p) e  ^  P^yttdt 
o 


/  >  f  ,  -(l-p)yt.. 

yp (1-p)  J  te  dt 

o 


P 

y (i-p)  * 


(30) 


Also,  the  probability  of  waiting  no  more  than  time  t  can  be 
obtained  as 


P(T<t)  =  1  -  P  (T>t )  =  1 


r  \  -d-p)yx, 

J  yp (1-p) e  K  dx 

t 


=  1 


pe-d-p)iit 


(31) 


If  instead  of  the  waiting  time  in  the  line  we  are  interested  in 


\ 


. 

' 
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the  total  length  of  time  spent  in  the  system  we  can  obtain  in  an 
analogous  manner  the  density  function  of  the  distribution  of  the  total 
time  in  the  system  of  the  (n+l)st  customer  as 


y  (yt) 


n  -yt 
e 


n! 


(32) 


which  implies  that  the  density  function  of  the  distribution  of  the  total 
time  in  the  system  of  any  customer  is 


U (1-p) e 


(1-p) yt 


(33) 


t  is 


Thus  the  probability  of  being  in  the  system  no  more  than  time 


1  -  e'(1'p)Mt 


and  the  expected  time  in  the  system  is 


(34) 


/  y ( 1-p ) e  U  P)uttdt 
o 


1 

yd-p) 


(35) 


2.2  The  Case  of  an  Infinite  Number  of  Channels 

Although  this  simple  model  is  useful  in  bringing  to  light  the 
types  of  results  which  can  be  obtained  it  is  not  very  realistic  from 
the  point  of  view  of  a  telephone  exchange  problem.  Let  us  modify  it  by 
assuming  that  instead  of  one  channel  we  have  infinitly  many.  In  this 


. 


■ 


21 


case  the  probability  of  a  call  being  initiated  during  a  small  interval  of 
time  At  will  still  be  AAt  +  o (At) ,  completely  independent  of  the 
number  of  occupied  channels,  whereas  the  probability  of  a  call  being 
terminated  during  a  small  interval  of  time  At  will  be  nyAt  +  o(At) 
where  n  is  the  number  of  occupied  channels . 

Thus  if  we  let  P  (t)  be  the  probability  of  n  occupied 

n 

channels  at  time  t  we  obtain  the  system  of  equations 

dP  (t) 

— —  =  -  ( A+nu) P  (t)  +  AP  (t)  +  (n+l)uP  (t)  , 
dt  n  n-1  n+1 

n  _>  1  / 

dP  (t) 

— §r —  =  -  AP  (t)  +  liP  (t)  .  (36) 

dt  o  1 

Proceeding  as  before,  we  again  wish  to  determine  the  generating 

function 


P  =  P(z,t)  =  I  P  (t)z  '  |z|  iL  1  >  (37) 

u  n 
n=o 


in  order  to  calculate  the  coefficients  of  zn  .  However,  in  this  case 
it  is  somewhat  easier  as  we  do  not  need  the  Laplace  transform. 

On  multiplying  the  first  equation  in  (36)  by  zU  and 
summing  over  the  appropriate  ranges  of  n,  including  the  second  equation. 


we  obtain 


- 


i 


.  . 


(1-z)  (- AP+y  — ) 

oZ 


(38) 


or ,  in  standard  form  , 


8P 

3t 


/  -j  .  3P 
(1-z)y  3^  = 


A(l-z)P  . 


(39) 


The  solution,  using  ordinary  methods,  is  obtained  as  follows: 
equation  (39)  has  the  associated  Lagrange  equations 


dt  _  dz  _  dP 
1  -  ( 1-z)  y  -A  ( 1-z)  P 


(40) 


Now 


dt 

~T 


dz 

-  (l-z)y 


yields 


u  (z  ,  t ,P) 


(1-z) 


c 


1 


and 


dz 

- (1-z) y 


dP 

-A  (1-z) P 


yields 


v (z , t ,P)  =  Pe 


- (A/y) z 


=  c. 


Hence  the  general  solution  is  given  by 


e(X/lJ)zg[(l-z)e-|Jt] 


P 


(41) 


. 
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If,  at  t  =  0,  i  channels  are  busy  we  have  that  P(z,0)  =  z1 


i  ( A/y) z  r  . ,  .  . , 

=>  z  =  e  g  [  (1-z)  e°  ] 


.  r  / -i  >  -Pt!  _  ~  (A/y)  (1-  (1-z)  e  Pt)  ,  ~yt  i 

:>  g[(l-z)e  ]  =  e  (l-(l-z)e  )  .  (42) 


When  this  is  substituted  in  (41)  we  obtain 


P(z,t)  =  e(X/y)Ze  (X/u)  (1  (1  z)e  5  (1-  (1-z)  e  yt) 1 


=  {exp  [-(A/y) (1-z) (1-e  yt) ] } [ 1- ( 1-z) e  Pt] 1  .  (43) 


Thus  we  have  obtained  an  expression  for  P(z,t)  and  therefore 
can  compute  P  (t)  be  calculating  the  coefficient  of  zU . 


However,  as  this  coefficient  is 


n 


p  <t>  =4 3  P<zit> 

n  ! 


ni  «  n 

9z 


z=0  ' 


(44) 


we  obtain 


P  (t) 
n 


-(A/y) (1-e  yt) 
nT~ 


l 


k=o 


n 

k 


) ( A/y) n  k 


i! 

(i-k) ! 


-yt,n+i-2k  -kyt 
(1-e  )  e 


for  n  <  i 


(45) 


and 


? 


- 
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P  (t) 
n 


-  (A/y)  (1-e  Ut) 


n! 


I 

k=o 


n  w,  ,  ,  n-k  i  ! 

k  > (Vw)  iwiT 


..  -iit.n+i-2k  -kyt 
(1-e  )  e 


for  n  >  i 


(46) 


by  applying  Leibnitz's  differentiation  formula 


dnf(z)g(z)  v  /  n  ,  ^(n-k)  (k)  ,  x 

- -  =  l  (  ,  )f  (z)g  (z) 


dz 


n 


(47) 


k=o 


to  P(z,t)  with  f  (z)  =  e^^^  ^  S  ^ Z  and  g(z)  =  [l-(l-z)e  ^Jt]  ^  . 


In  order  to  find  the  average  number  of  channels  occupied  at 
time  t  ,  multiply  the  first  equation  of  (36)  by  n  and  add  over  the 
appropriate  ranges  of  n  .  This  gives 


y  nP'  (t)  =  M'  (t)  =  A  -  yM (t) .  (48) 

i  n 

n=l 


If  we  assume  at  t  =  0,  i  channels  are  busy,  then  M(0)  =  i 

and 


M  ( t) 


(A/y ) (1-e  Pt)  +  ie 


(49) 


It  is  interesting  to  note  that  in  the  special  case  i  =  0  we 


have 


? 
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P  (t) 
n 


-(A/y)(l-e  yt)  ,w  .  n  -yt  n 
e  _ (A/y)  ( 1-e  ) 

n! 


(50) 


that  is,  the  P  (t)  are  0iven  exactly  by  the  Poisson  distributions  with 
mean  M(t)  . 

To  determine  the  equations  for  the  probabilities  which 

describe  the  steady  state,  as  before  we  set  the  time  derivatives 
equal  to  zero  and  eliminate  time  from  the  system  of  equations  (36) , 
obtaining 


/ 


(A+ny)p  =  Ap  +  (n+l)yp  ,  n  >  1  ,  (51) 

n  n-1  n+1  — 


which  gives  by  induction 


P 


n 


=  Po(A/y)n/n! 


Pn  = 


e  (X/y)(A/y)n 

n! 


(52) 


Thus  the  limiting  distribution  is  a  Poisson  distribution  with  parameter 
A/y  .  It  is  independent  of  the  initial  state. 


. 


7 
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2.3  The  Case  of  a  Finite  Number  of  Channels — Delay  System 

In  the  above  example  we  are  not  concerned  with  waiting  lines 
as  we  have  an  infinite  number  of  channels.  However,  the  situation  is 
still  somewhat  unrealistic.  If  we  assume  instead  that  we  have  a 
finite  number  of  channels,  say  c  ,  and  that  if  all  channels  are  busy  a 
single  waiting  line  forms  and  calls  enter  the  channels  when  they  are 
vacant  on  a  f irst-come-f irst-served  basis ,  we  will  approximate  the 
situation  of  a  long  distance  switchboard. 

Again  by  letting  P  (t)  be  probability  of  n  calls  in 

the  system  at  time  t  we  obtain  the  system  of  equations 


dP  (t) 
o 

dt 


-  AP  (t)  +  yP.  (t)  , 
o  K  1 


dP  (t) 
n 

dt 


(A+ny)P  (t)  +  AP  ■,  (t)  +  (n+l)yP  (t), 
n  n-1  n+I 


(53) 


1  <  n  <  c  , 


dP  (t) 
n 

dt 


(A+cy)P  (t)  +  AP  ,  (t)  +  cyP  (t), 
H  n  n-1  n+1 


c  <  n  . 


We  begin  the  solution  of  these  equations  in  much  the  same 
manner  as  those  in  Erlang's  model  by  considering  the  generating 
function 


=  P (z,t)  =  J  P  (t) z  , 
u  n 
n=o 


P 


z 


<  1  , 


(54) 


'  •  :  I 
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and  taking  the  Laplace  transform. 


Rewriting  the  system  (53)  in  the  form 


P'(t)  =  -  (X+cy ) P  (t)  +  cyp  (t)  +  cyP  (t)  -  (c-l)yp  (t), 
o  o  o  1  1 

P'  (t)  =  -  ( X+cy) P  (t)  +  (c-n)yp  (t)  +  XP  ft) 
n  n  n  n-1 

+  cyP  (t)  -  (c-n-l)yP  (t)  ,  1  <  n  <  c  , 

n+1  n+1  — 

P'  (t)  =  -  (X+cy) P  (t)  +  XP  (t)  +  cyP  (t)  ,  c  <  n, 

n  n  n-1  n+1  — 

multiplying  each  equation  of  (55)  by  z*1  and  summing  gives 


9P 

3t 


P-P  (t) 

(X+cy)P  +  XzP  +  cy  - - -  +  cyP  (t) 

z  o 

(c-1) y (1-z) P1 (t)  -  (c-2) yz (1-z) P2 (t)  -  ... 


(c-n)yzn  1 (l-z)P  (t)  - 

n 


-  yz°  2  (l-z)P  (t) 
c-1 


Taking  the  Laplace  transform  of  this  expression  and  solving 
for  P*(z,s)  gives 


P* (z,s) 


c-1 

z^+^~  (1-z) y  1  (c-n)znp*(s) 
n=o  n 


if  we  assume  P. (0)  =  1  . 

l 


(55) 


(56) 


sz- (1-z)  (cy-Xz) 


(57) 


. 

■ 
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Since  P  (z,s)  converges  inside  and  on  the  unit  circle  for 
Re(s)  >  0,  the  zeros  of  the  numerator  must  correspond  with  those  of  the 
denominator  in  this  domain. 

The  zeros  of  the  denominator  are 


a.  (s)  = 
k 


2  b 

A+cy+s± [ ( A+cy+s)  -4Acy] 
2A 


k  =  1,2 


(58) 


where  a^(s)  has  the  positive  sign  before  the  radical. 

Again,  applying  Rouche's  Theorem  to  the  denominator 

sz  -  (1-z) (cy-Az) ,  we  find  that  a  (s)  is  its  only  zero  inside  the 

•  c-1 

i+1  v  n  * 

unit  circle  and  thus  the  numerator  z  -  y(l-z)  l  (c-n)z  P  (s) 

n=o 

must  vanish  for  z  =  ,  that  is 


l  (c-n)a  P  (s) 
L  2  n 

n=o 


i+1 

_^2 _ 

U(l-a2) 


(59) 


Now  in  order  to  obtain  P  (s) ,  and  thus  P  (t) ,  for 

n  n 

0<n<c-l  we  require  c  equations  involving  the  c  unknowns 

p* (s) ,  0  <  n  <  c  -  1  .  These  c  equations  are  (59)  and  the  first 

n  —  — 

c  -  1  equations  of  (53)  which  are 

P'  (t)  =  -  AP  (t)  +  yP  (t)  , 
o  o  1 

P '  ( t)  =  -  ( A+ny)  P  (t)  +  AP  (t)  +  (n+1)  yP  (t)  ,  (60) 

n  n  n-1  n+1 


1  <  n  <  c 


2. 


. 


cx  -P 
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To  solve  these  equations  we  introduce  the  generating  function 


c-2 

Q(z,t)  =  l  Pn(t)zn  .  (61) 

n=o 


Multiplying  the  second  equation  of 


(60) 


by 


n 


z 


and  summing 


gives 


•XS  -  UZ  ||  +  y  |2  +  Azg  -  AzC_1Po_2  + 


(62) 


9Q 

3t 


=>  -  qr-  +  y  (i-z) 


30  ,  ,  _  .  .  c-1 

=  A(l-z)Q  +  Az  P 
3z 


c-2 


.  ,  .  c-2 

-  y  (c-1)  z  P 


c-1 


(63) 


dt  _  dz  _  _ _ dQ _ 

-1  y(l-z)  A  (1-z)  Q+Azc--*-P  -y  (c-1)  zc_^P 

c-2  c-1 


(64) 


The  first  equation  gives 


c1eyt  =  (1-z) 


C-l  =  (l-z)e 


-yt 


=>  z  =  1  - 


Cl6 


yt 


(65) 


Using  this  value  of  z  in  the  third  member  and  taking  the 
equation  in  the  first  and  third  members  gives 


M  +  xc/h?  -  u (c-i)  -  mi-c/v\2 


(66) 


with  the  solution 
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Q  =  e 


- ( A/y) c1eyt  t 


r  r  r  ,  , .  yxv  c-2  .  .  .  ,,  yx%  c-r  .  .  . 

{{  [y (c-1)  (1-c  e  )  P  (x)-A(l-c  e  )  P  (x) ] 

J  n—  I  n—  ) 


yx, c-1 


o 


1  c-1 1  '  c-2 


(A/y)CleUX 
e  dx+c^} 


(67) 


In  order  to  obtain  c^  =  f(c^)  we  must  consider  two  cases. 
First  of  all  if  the  initial  number  in  the  system  i  <_  c  -  2, 

Q( z,0)  =  z1 


.  jr  / 1  ,  i  (A/y)  (1-z) 

=>  f(l-z)  =  z  e 


(68) 


as  f(cx)  =  f  [  ( 1-z)  e  yt] 


f [ ( 1-z) e  yt]  =  [1- (1-z) e  yt]X  exp  [(A/y)(l-z)e  Ut]  (69) 


==>  Q  ( z ,  t)  =  e  (A/u)(1  z)  /  {y(c-l)  [l-(l-z)e  y(t  x)  ]  C  ^^(x) 


o 


-  A [1- (1-z) e  y(t  X) ]C  XP  (x) } 

c-2 


r , ,  ,  ,  ,,  x  -y  (t-x) 
exp  [ (A/y)  (1-z) e  ] dx 


+  exp  [-(A/y)  (1-z)  (1-e  y ) ]  [ 1- ( 1-z) e  yt] 1  .  (70) 


Secondly,  if  the  initial  number  in  the  system  i  >  c  -  2, 
Q(z,0)  =0  and  hence  c^  =  0  . 


Continuing  the  problem  for  i  >  c  -  2 ,  we  use  the  result; 


3Erdelyi,  A.  (ea.) :  "Tables  of  Integral  Transforms",  vol.  I, 
McGraw  Hill  Book  Co.,  Inc.,  New  York,  1954,  page  147,  (40). 


. 
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—  IJT/j 

f  ~sw  / 1  v  -yw  c-2  (A/y)  (1-z)  e 
J  e  [1- (1-z) e  ]  e  dw 

o 

00  - TT 

=  J  e~(s/y)Y[l-(l-z)e_Y]C"2e(Vy)(1'2)e  & 

=  y  r((s/y)+l)  ^l[y  (c_2)  ,:L'  (1_z)  '  (A/u)  (1"z)  ] 


where 


(J>1  ( a,  3  ,  Y  ,x  ,y) 


00  oo 

(a)m+n(g)n  m  n 
(y)m+nm!n!  X  Y  ' 

n=o  m=o 


l  l 


(a)  =  1  and  (a)  =  a (a+1) . . .  (a+n-1)  . 

o  n 


This  gives 


Q* (z  ,s) 


e 


(A/y) (1-z) 


{ (c-l)P*  . (s) 
c-1 


T (s/y) 

T ( (s/y) +1) 


eft  [s/y  (c-2)  ,1,  (1-z)  ,  (A/y)  ( 1-z)  ]  -  ( A/y )  P*  _(s) 
1  c-2 


T (s/y) 

K  (s/y)+l) 


(j>1  [s/ y ,-  (c-1)  fl,(l-z)  ,  (A/y)  (1-z)  ]  } 


(71) 


as  the  first  expression  of  (70)  is  a  convolution. 


We  have 


P 


* 

c-2 


(s) 


i  8°-y 

(c-2) !  c-2 

dZ 


z=0 


(72) 


? 


) 


i'v\  >:  I  .  ;  •  •  X.'C.l  '  -  '  .C  3~7  ±  'Jf.r  TV. 
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c-2 

*  .  r  *  .  ,  n 

as  Q  (z,s)  =  )  P  (s)z 

u  n 
n=o 

Letting  f(z)  =  e  and  g(z)  =  <J)  [s/y (c-2)  ,  1 ,  (1-z)  , 

(A/y)  (1-z)  ]  ,  applying  the  Leibnitz's  differentiation  theorem  gives 


1 _  3°  2f(z)g(z) 


(c-2) !  c-2 

3z 


z=0 


c-2 

1  v  ,c-2 


l  (CVZ)  (A/y)C  k  2e  X/U  cj)1  ,-  (c-2)  ,l,(l-z)  ,(A/y)  (1-z)]  |_ 


(c-2)  !  L  k 
k=o 


dz 


k  1  p 


z=0 


(73) 


Thus 


p*  (S)  =  _ L(s/y) _ 

c-2V  ;  (c-2) !f( (s/y)+l) 


c-2 

,  *  v  c— 2 

{(c-l)Po_1(s)  l  (  k  ) 

k=o 


( A/y)  C  k  2e  X/y  <f>1  [s/y ,-  (c-2)  ,1,  (1-z)  ,  (A/y)  (1-z)  ]  |  z=Q 

dz 


c-2 


(A/ti)p*-2(s)  l  t/xVw) 
k=o 


c-k-2  -A/y 
e 


dz 


k  T1 


<K  [s/y , -  (c-1)  ,1,  (1-z)  ,  (A/y)  (1-z)  ]  lz=Q> 


(74) 


Solving  for  pc_2^s^  yields 


I 
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P 


* 

c-2 


(s) 


r  r  (s/yi)  (c-1) 

(c-2)  !  T  (  (s/y)  +1) 


PC-1(S)  l  (°k2) (x/y) 

k=o 


c-k-2  -A/y 
e 


<f)  [s/y,- (c-2)  ,1,  (1-z)  ,  (A/y)  (1-z)  ]  |  }  * 

,  K  1  Z=U 

dz 


(1  +  (A/y) 


T  (s/y) 

(c-2)  !  T  (  (s/y)  +1) 


c-2 

l 


k=o 


(C,  2)  (A/y) 
k 


c-k-2  -A/y 
e 


<#)  [s/y,- (c-1)  ,1,  (1-z)  ,  (A/y)  (1-z)  ]  |  _}  • 

n  K  i  Z=U 

dz 


We  also  have  in  exactly  the  same  manner  as  above  for  n  < 


P*(s) 


T  (s/y) 


n!  T  (  (s/y)  +1) 


{  (c-l)P  .  (s) 
c-1 


n 

l 

k=o 


n 

k 


n-k  -A/y 
)  (A/y)  e 


<J>,  [s/y,  -  (c-2)  ,1,  (1-z)  ,  (A/y)  (1-z)  ] 

dzk 


,  \  v  /  n  \  ,w  xn-k  -A/y 

-  ( A/y)  P  (s)  ),  (  v  )  (A/y)  e 

C”  A  .  K 

k=o 

X 

— —  cj)  [s/y,  -  (c-1)  ,1,  (1-z)  ,  (A/y)  (1-z)  ]  |  }  . 

k  1  z=0 

dz 


Thus  as 


\  (c-n)a  P  (s) 
L  2  n 

n=o 


i+1 

_^2 _ 

y  (1-a2) 


(75) 


-  2 


(76) 


(77) 


we  have 


>  -('■ 


\ 


34 


a. 


i+1 


P  ,  (s) 
c-1 


U (l-a2) 


c-2 

*  +  1  (c-n) 

n=o 


* 

P  (s) 


n  n 


a. 


2  * 
c-1 


(78) 


P"  , (s) 


and  thus  we  have  our  c  equations  involving  the  c  unknowns  , 

0  <  n  <  c  -  1. 


In  order  to  determine  P  (s)  ,  n  _>  c,  we  will  again  consider 


the  generating  function 


i+1 


c-1 


P  (z,s)  = 


n=o 


-y(l-z)  £  (c-n)znp*(s) 

u  n 


sz  -  (1-z)  (cy-Az) 


(79) 


and  determine  the  coefficients  of  zU .  To  do  this  write  the  denominator 
in  the  form  -A(z  -  a^)  (z  -a  )  and  apply  Leibnitz's  theorem  to  the 
product  (z-a  )  ^(z-a^)  ^ .  We  have 


1  dm  ,  -1.  -li 

— - (z-a.)  (z-a  ) 

m!  m  1  2  z=0 

dz 


1 


m 


.1(7)  (-D  k!  (-a  ) 
m!  ^  k  ' 

k=o 


- (k+1)  .m-k 
y  (-D  (m-k)  !  (-a2) 


-  (m-k+1)  (80) 


m 


V  ,  v - (k+1)  - (m-k+1) 

L  (a±)  (a2} 

k=o 


m 


k=o 


V  -1  -m-1 .  .  . . 

I  a1  «2  (“/V 


X  i-fc/ay 


m+1 


m  1-  (a„/an ) 
cya2  2  1 


(81) 


'  . 


? 


. 


35 


Now 


and  thus 


p  (s)  =  — - - P  (z,s) 

n  n !  .  n  1 z=0 

dz 


p*(S)  =  4  1  ( i ) 

n  n!  .  u  k 

k=o 


,  k  /v  ,n-k  .  . 
d  A(z)  d  B(z) 


dz 


dz 


n-k  1 z=0 


c-1 


c-1 


where  A(z)  =  z^+^  -  y  £  (c-j)z"^P*(s)  +  y  \  (c-  j )  z3  +  3~'P*  (s)  and 

j=o  11  j=o  3 


,  s  -1,  -1  A+cy+s  cy 

B  (z)  =  (-1/A)  (z-a1)  (z-a^)  as  +  a2  =  - - -  ,  =  —  and 

s  =  -A  (l-a1)  (l-a2)  , 


:  >  P 


n 


c-1  ,n-i  /  x 

,  ,  -y  v  ,  n  .  .  .  x  .  ,„*  /  x  d  B (z)  | 

(S)  =  -7  )  .  )  (C-U  :!P.  S  - : n 

nl  ,L  i  l  n~ 3  2=0 

j=o  dz 


c-1  nn— ( j+1 ) ^ ,  \ 

(j;1MH»ij+i»^(s, g-.n-.(^iz=o 

j=o  dz 


+ 


n- ( i+1) 
d  B  (z) 


(n-i-1)!  ,  n-(i+l)  z=0 

dz 


c-1 

=  (y/A)  l  [(c-j)P*(s) 

j-o 


A  1_(a2/al) 


n- j+1 


n-j  1-a  / a 
cya2  2  1 


,  1-  ( a  /  a  ) n  D 

-  (c-:)P.(s)  - — - - - - ] 

3  _ n-j-1  l-a0/an 


cya„ 


2  1 


1  l-(a2/a1) 


n-i 


n-i-1  l-a_/a-. 

cya  2.  1 


(82) 


(83) 


(84) 


for  n  >  c  . 


(85) 


.'fS 
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Thus  we  have  obtained  P*(s)  for  n  =  0,1,2,...  and  now  may 
obtain  P  (t)  for  n  =  0,1,2,...  by  considering  the  inverse  transform. 

In  the  same  way  as  before  we  obtain  the  equations  describing 
the  steady  state  as 


ApQ  =  up  , 


(A+ny)p^  =  Ap^_1  +  (n+l)ypri+1  >  1  £  n  <  c  , 


n-1 


n+1 


(86) 


(A+cy)pn  =  Apn_]_  +  cypn+1  , 


n  >  c  . 


By  induction  we  find  that 


(A/y) 

p  =  P  - — 

n  o  n ! 


n 


n  <  c 


P  =  P 
n  o  ,  n-c 
c !  c 


(A/y) n  PQ<Vcy)ncC 


c  1 


n  >  c  • 


(87) 


r  n 

However,  the  series  i  —  converges  only  if  A/y  <  c  . 

^o 

Hence  a  limiting  distribution  cannot  exist  when  A  >_  cy  .  In  this  case 
p^  =  0  for  all  n  ,  which  means  that  gradually  the  waiting  line  grows 
over  all  bounds. 


On  the  other  hand,  if  A/y  <  c,  using  \  p^  =  1,  we  obtain 


n=o 


c-1  ,,  .  ,  n  00  n  c 

i  =  P  i  l  +  l  e-f-i 

o  L  nl  L  c! 

n=o  n=c 


(88) 


- 
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where  p  =  —  , 

cy 


= 


0-1  /  sn 

V  (cp)  +  (cp) 


n=o 


n!  c ! ( 1—  p ) 


(89) 


The  expected  number  in  the  system  is  given  by 


°°  c  /  -,  /  , n  00  ,  .  n  n 

L  =  l  np  =  p  y  +  p  y  n  <X/cU)  c 

L  n  *o  L  n!  o  L  n  c! 

n=o  n=o  n=c+l 


(90) 


p(cp) 

=  cp  +  - -  p 


c! (1-p) 


2  o 


(91) 


whereas  the  expected  number  in  the  queue  is 


L 

q 


CO  Q 

y  /  x  p (cp) 

I  (n-c)p  =  - 

n=c+l  n  c ! (1-p ) 


2  Po 


(92) 


Again  in  this  example  the  situation  may  occur  where  an 
arrival  will  have  to  wait.  The  probability  of  this  occuring  is  clearly 
the  probability  that  c  or  more  channels  are  occupied  which  is 


I 

n=c 


n  c 
(p )  c 


c 

pc  c 

12 _  .  P 

c !  1-p 


p  (cp)C 
o _ 

c !  (1-p ) 


(93) 


To  determine  the  probability  P(T>t)  of  waiting  longer  than 
time  t  consider  the  probability  of  the  nth  customer  in  the  line 
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waiting  longer  than  time  t.  This  is 


P{T +T  +.  .  .+T  >  t> 

12  n 


(94) 


where  T_^  is  a  random  variable  giving  the  amount  of  time  all  c 
channels  are  occupied  before  any  of  them  becomes  free. 


Thus  p{t^  >  t}  =  the  probability  that  none  of  the  c 

—  cut 

become  free  in  time  t  which  is  e  and  therefore  T.  has 

x 

exponential  distribution  with  mean  cy  . 


channels 

an 


is 


cy 


cy-ic}) 


The  characteristic  function  of  this  exponential  distribution 
and  thus  the  characteristic  function  of  the  distribution  of 


the  waiting  time  of  the  nth  customer  in  the  line  is  (— C-^-.-,-) n  ,  thus 

cy-i4> 


giving  a  density  of 


n-1  -cyt 
cy(cyt)  e 

(n-1) ! 


(95) 


which  implies 


W  (t)  =  the  probability  that  the  nth  customer 
n 

in  the  line  is  served  by  time  t 


t  n-1  -cyx 

I  (cy) (cyx)  e  dx 


o 


(n-1)  ! 


(96) 


The  density  function  of  the  distribution  of  the  waiting  time  of 


any  customer  in  the  line  is  given  by 


O 

■ 


J 


.  ;'b 


. 


39 


P 


n=l 


c+ (n-1) 


,  .  .  n-1  -cyt 

(cy) (cyt)  e 

(n-1)  ! 


p  c  00  ,  /  \  t  ^ n-1  -cyt 

o  r  .  c+ (n-1)  (cy) (cyt)  e 
-  L  (P)  - 


c! 


n=l 


c.  ,  -cyt  c 
p  c  (cy) e  p 
o  r  (pcyt) 


(n-1)  ! 

n-1 


C  ! 


y  ipcyt; 

\  (n-1) ! 


n=l 


p  cC  (cy) p° 
o _  -cyt(l-p) 

c!  8 


(97) 


and  thus  we  can  calculate  the  probability  that  an  arriving  customer  will 
wait  longer  than  time  t  as 


P (T>t )  = 


p  c  (cy)  p  °o 
o  f  -cyx ( 1-p )  „ 

- J  e  dx 

c !  J 


PoC  P  -cyt (1-p) 


c!  (1-p) 


=  e-^td-P)p(T>0) 


(98) 


Approaching  this  problem  from  a  different  point  of  view,  we 
clearly  have  W^(t)  equal  to  the  probability  that  any  of  the  c  channels 
completes  service  in  time  t  and  thus 


t 

W  (t)  =  /  W  (t-x)d  W n  (x) 
n  J  n- 1  x  1 

o 


=  1 


-cyt 

e 


I 


i=o 


(cyt) 1 

i! 


(99) 


by  induction. 
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The  probability  of  waiting  less  than  time  t  is  thus 


P(<  t) 


c-1 


n=o 


P  +  L  P  ,  /  ,  (t) 

n  u _  c+(n-l)  n 

n=l 


(100) 


where  W  (t)  is  given  above, 
n 

2.4  The  Case  of  a  Finite  Number  of  Channels — Loss  System 

In  concluding  this  first  part  we  alter  the  delay  system  just 
discussed  by  assuming  that  a  call  is  completely  lost  to  the  system  if 
all  c  channels  are  busy  (that  is,  no  waiting  line  is  formed).  This  is  the 
situation  of  the  trunks  that  are  available  to  a  long-distance  operator 
when  she  places  a  call  for  a  customer  or  those  available  to  a  certain 
group  of  subscribers  for  the  use  of  their  telephones. 

Again  by  letting  P  (t)  be  ^he  probability  of  n  calls  in 
the  system  at  time  t,  we  obtain  the  system  of  equations 

P'  (t)  =  -  IP  (t)  +  yP.  (t)  , 
o  o  1 

P'(t)  =  -  (A+ny ) P  (t)  +  AP  (t)  +  (n+l)yP  (t)  , 
n  n  n-1  n+1 

1  <  n  <  c  -  1  , 


P'  (t)  =  -  cyP  (t)  +  AP  (t)  .  (101) 

c  c  c  —  1 


These  differ  from  the  previous  situation  only  in  the  last 


equation. 
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The  equations  describing  the  steady  state  are 


Apo  =  PPl  , 


(A+ny)p^  =  ^Pn_1  +  (n+1)lJPn+1  '  1  <_  n  <_  c  -  1  , 


c|jpo  =  Apc-1  ' 


and  as  before  we  have 


p©(  Vy) 

Pn  n ! 


n 


Since  Y  p  =  1, 
u  n 


n=o 


which  implies 


£  1AM^  =  q  , 

o  L  n! 
n=o 


Po  = 


1  +  A/y  +  (A/y)  /2!  +  ...  +  (A/y)  /cl 


and  thus 


Pn  = 


Pn/ni 


1+  p  +  p  /2!  +  ...  +p  /cl 


(102) 


(103) 


(104) 


(105) 


(106) 


where  p  =  A/y  . 
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Since  the  probability  that  a  call  is  lost  equals  the 
probability  that  all  c  channels  are  busy  we  have  the  probability  of 
a  lost  call  as 


pC/ci 


1  +  p  +  p  /2!  + 


+  pC/c ! 


(107) 


which  is  called  the  Erlang  loss  formula. 

As  c  ->  00  we  have  the  case  of  an  infinite  number  of  channels. 
Letting  c  ->  00  in  (106)  we  obtain  the  Poisson  distribution 


Pn  = 


Pn/n! 


2  n 
l+p+p  /2 !  +  ...  +p  /n!  + 


n 

PUe'P  , 

n ! 


p  =  X/y  , 


(108) 


which  agrees  with  out  previous  result. 


As  clearly  there  are  no  waiting  times  in  this  situation,  we 


have  completed  the  first  part  of  the  thesis  and  are  ready  to  begin  the 
analysis  of  the  telephone  usage  of  a  group  of  local  subscribers. 


i 


CHAPTER  III 


AN  ACTUAL  TELEPHONE  EXCHANGE  SITUATION 


3.1  A  Description  of  the  Data 

The  data  that  have  been  used  in  this  analysis  were  recorded  by 
a  pen  recorder  chart  on  the  calls  terminating  at  the  699  77 —  Sherwood 
Park  telephones  from  9  A.M.  to  5  P.M.  over  a  period  of  about  eight 
days.  The  actual  days  and  times  considered  are  as  follows: 


Thursday , 

October 

2nd , 

1969 

9:00 

A.M. 

to 

5:00 

P.M. 

Friday , 

October 

3rd , 

1969 

9:00 

A.M. 

to 

5:00 

P.M. 

Monday , 

October 

6th, 

1969 

9:00 

A.M. 

to 

5:00 

P.M. 

Tuesday , 

October 

7  th, 

1969 

n  9:00 
and  „ , 

11:00 

A.M. 

to 

10:39 

A.M. 

A.M. 

to 

5:00 

P.M. 

Wednesday , 

October 

8  th, 

1969 

9:00 

A.M. 

to 

5:00 

P.M. 

Thursday , 

October 

9  th, 

1969 

9:00 

A.M. 

to 

5:00 

P.M. 

Friday , 

October 

10th, 

1969 

9:00 
and  2:30 

A.M. 

to 

12:00 

P.M. 

P.M. 

to 

5:00 

P.M. 

Tuesday , 

October 

14th, 

1969 

9:00 

A.M. 

to 

1:00 

P.M. 

Since  the  699  77 —  telephones  are  all  business  phones  which 
operate  on  a  rotary  hunting  principle,  the  times  and  days  considered 
give  a  reasonably  accurate  description  of  their  use. 


When  one  dials  a  number  beginning  with  699  77 —  it  is 
necessary  to  make  a  connection  through  one  of  the  ten  trunks  represented 
by  the  first  ten  lines  (lines  1  to  10)  of  the  pen  recorded  chart.  (See 
the  example  on  page  44) .  At  the  instant  when  the  fifth  digit  in  the 
number  has  been  dialed,  a  connection  is  made  with  one  of  the  ten  trunks 
and  the  pen  corresponding  to  that  trunk  will  jump  upward,  remaining  in 
its  new  position  until  the  completion  of  the  call.  A  call  that  is 
attempted  when  all  channels  are  busy  is  lost.  Clearly,  we  have  the 
fourth  situation  discussed  in  Chapter  II  -  that  of  a  finite  number  of 
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Sample  of  the  Pen  Recorder  Chart 


1  °  '• 


45 


channels  available  to  a  group  of  subscribers  where  a  call  is  lost  if  all 
channels  are  occupied. 

Thus,  by  again  assuming  that  all  calls  arrive  independently 
of  each  other  and  of  the  number  in  the  system  and  that  the  length  of  each 
call  is  independent  of  every  other  call,  we  can  compare  the  distributions 
calculated  from  the  interarrival  times  and  the  busy  times  of  the  given 
data  with  the  exponential  situation  discussed  in  Chapter  II. 

The  second  ten  lines  (lines  11  to  20)  on  the  chart  show  both 
the  originating  and  terminating  calls  on  particular  Sherwood  Park 
telephones  and  have  not  been  considered  in  this  analysis. 

3.2  The  Method  of  Analysis 

The  data  were  considered,  first  of  all,  on  an  hourly  basis  and 
the  number  of  interarrival  and  busy  times  of  the  various  lengths  were 
recorded.  Since  it  was  not  feasible  to  measure  lengths  on  the  chart  more 
accurately,  they  were  grouped  into  intervals  of  length  ten  seconds,  which 
is  represented  by  the  smallest  division  on  the  recording.  These  results 
are  given  in  Tables  1.1  to  1.16  of  Appendix  I. 

As  the  object  of  this  analysis  is  to  see  if  these  distributions 

are  actually  exponential,  it  was  first  necessary  to  estimate  the  parameters 

X  and  u  where  —  and  —  are  the  average  interarrival  and  busy  times 

X  U 

respectively . 

Thus ,  if  T^  is  the  interarrival  time  in  seconds ,  we  assumed 
that  T  has  an  exponential  distribution  with  parameter  X  where  X 
is  unknown.  Letting  xg ,X2 ' * ' * ,xn  be  N  indePendent  observations  of 


■ 


. 
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T  with  joint  density  function 


N 

-A  £  x . 

,N  i=l  1  N  -ANx 

A  e  =  A  e 


we  have  the  maximum  likelihood  equation 


{N  In  A-A  Nx}  =  0  , 
dA 


which  has  the  unique  solution 


A  = 


x 


From  the  fact  that  the  second  derivative  of  the  log- likelihood  function 


-  {N  In  A-A  Nx}  =  — 
dA2  A2 


is  everywhere  negative,  we  know  that  the  above  solution  is  the  maximum 
likelihood  estimate. 


Now,  if  n  is  the  number  of  observations  in  [r,r+l)h, 

r  s  , 

r*  _  O 

n  ,  and  S  =  )  rn  ,  we  have  X,  =  ( —  +  — ) 

r  o  u  r  1  N  2 

r>o  r>o  o 


=0,1,2,...,  N  =  y 

o  u 

»  T, 


where  =  — —  and  h  -  the  unit  of  grouping  =  ten  seconds.  Thus,  for 

1  h 

this  grouped  data. 


T  =  (—  +  -q-)  *h 
1  N  2 
o 


giving 


< 


•  fit 
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{<r  +  !>-h} 

o 


However , 


P{rh<T1< (r+1) h}  =  0^ (1-0)  ,  r  =  0,1,2,...  , 


where  0^  = 
A  we  have 


-Ah 


e 


r 


and  on  substituting  the  maximum  likelihood  value  of 


A  -Ah  o  1-1 

®1  =  6  =exP-<iT+I> 

o 


The  situation  was  not  quite  so  straightforward  when  dealing 
with  the  busy  times  .  The  problem  arises  as  the  busy  time  shown 

by  the  pen  recorder  is  actually  the  sum  of  two  random  variables,  X 
and  Y.  X  gives  the  time  required  to  dial  two  digits  (approximately 
five  seconds),  the  ringing  time  (approximately  nine  seconds,  if  we 
assume  two  rings)  and  time  necessary  to  identify  oneself  and  call  the 
desired  party  to  the  telephone  and  Y  gives  the  actual  conversation 
time.  We  have  assumed  that  Y  is  exponential  but  we  have  no 
information  concerning  the  distribution  of  X  and  thus  have  assumed  it 
constant  and  equal  to  thirty  seconds.  This  choice  is  completely 
arbitrary.  The  actual  value  may  vary  greatly  from  call  to  call  and 
cause  considerable  error  in  the  calculations. 


Letting 


Y 


T  -ah 
2 


Y 

h 


r 


a 


h 


a  =  3, 


h  =  10, 


!<■'  •>  • 


-•  •> 


n 


■: )  >  1  >  ‘  •  'T 


t'ft  ■’ • 
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where  the  conditional  probability  density  function  of  Y  given  Y  >  0 
is 


ye  Uy,  y  _>  0  , 

0  ,  y  <  0  , 

we  have,  by  the  same  type  of  analysis  used  for  the  interarrival  times, 
that 


giving 


a 


{<n  i  _a)'h} 

a 


-1 


where  N  =  J  n  and  S  =  £ 


a  “  r 
r>a 


rn 

a  “  r 
r>a 


Also , 


p{rh<Y<  (r+1)  h}  =  r  = 


where  0^ 


-yh  .  . 
e  ,  giving 


~  -yh  ,  a  1  -1 

02  =  e  =  exp  -  (—  +  -  -a) 

a 


The  maximum  likelihood  estimates  of  A  and  y  for  each  hour 
of  each  day  are  given  in  Tables  1.17  and  1.18  of  Appendix  I.  The 
corresponding  estimates  of  0^  and  0^  are  given  in  Table  1.19  and 
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1.20  of  Appendix  I. 


After  this,  the  total  observed  frequencies  for  each  one  hour 
period  were  calculated  and  the  results  tabulated  in  Tables  2.1  and  2.2 
of  Appendix  II.  Again  the  maximum  likelihood  estimates  of  A  and  P 
and  the  corresponding  0^  and  0^  were  obtained.  These  results  are 
given  in  Tables  2.3  and  2.4  of  Appendix  II. 


Finally,  using  these  estimates  for  0^  and  0^,  x  tests  were 
done  on  the  hourly  results,  facilitated  by  the  use  of  the  computer. 

These  results  are  given  in  Table  2.5  of  Appendix  II. 


Here  the  x  values  were  obtained  for  the  interarrival  times 


by  comparing  the  observed  frequencies  n^  with  the  expected 

frequencies  0(1-0_,)N  ,  r  =  0,1,2,...,  where  N  =  >  n 

1  1  o  o  u  r 

r>o 


The 


calculations  were  pooled  for 


r  =  9  and  10  , 
r  =  11  and  12  , 
and  r  >  13  , 


giving  10  degrees  of  freedom. 


The  computing  of  the  X  values  for  the  Y  distribution  was 

again  more  difficult.  Here  the  observed  frequencies  n  ,  r  =  3,4,5,..., 

/\  ^ 

were  compared  with  the  expected  frequencies  09  (1-0O)N^  where 
r 1  =  r-3 ,  r 
pooled  for 


=  0,1,2, ... ,  and  N  =  \ 

r>3 


n 


2  ~2'"3 

The  calculations  were 


18 

and 

19, 

20 

and 

21, 

r 


T  T 


. 

■ 
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r  =  22 
r  =  24 
r  -  26 
r  =  28 
r  =  31 
r  =  34 
and  r  >  38 


and 

23/ 

and 

25, 

and 

27, 

29 

and 

30 

32 

and 

33 

35, 

36 

and 

giving  22  degrees  of  freedom. 


The  graphical  representations  of  Appendix  III  give  a  further 
comparison  of  the  observed  and  expected  frequencies. 


3.3  Results  and  Conclusions 

a)  Throughout  this  analysis  it  has  been  obvious  that  the  results 
obtained  for  the  interarrival  times  during  the  lunch  hour,  from  12  noon 
to  1  P.M.,  differ  substantially  from  the  others. 

Considering  the  maximum  likelihood  estimates  given  in 
Table  2.3  of  Appendix  II,  the  value  of  A  for  this  time  is  certainly 
smaller  than  the  others.  This,  of  course,  should  be  the  case  as  the 
number  of  arrivals  during  the  lunch  hour  is  lower  and  thus  the  expected 
interarrival  time  is  greater . 

b)  The  probability  of  a  call  being  lost  is  equal  to  the  probability 
that  all  of  the  ten  channels  are  occupied. 

Altogether,  57  hours  9  minutes  =  3,429  minutes  =  20,574 
"ten  second  intervals"  were  analysed.  Of  this  time,  all  channels  were 
occupied  simultaneously  for  80  "ten  second  intervals". 
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Thus,  on  the  basis  of  this  data,  the  probability  of  a  call 
being  lost  is 


80 

20,574 


.0039  . 


c)  In  the  second  chapter  we  discussed  the  quantity  —  ,  the  traffic 

I* 

intensity.  In  order  to  obtain  an  approximate  value  of  —  for  our 

h* 

A  /S 

data,  average  values  of  X  and  y  have  been  calculated  from  Table  2.3 
of  Appendix  II,  giving 


X  =  .026,  y  =  .0082 

A 

and  —  =  3.15  • 

A 

y 

Since  we  are  dealing  with  ten  channels ,  the  expected  number 
of  arrivals  at  each  channel  will  be  ~  of  the  total  expected  number, 
implying  that  the  expected  interarrival  time  at  each  channel  will  be 

ten  times  greater,  that  is,  10  x  -  .  Thus  the  traffic  intensity,  or  the 

/\ 

X 

ratio  of  expected  busy  time  to  expected  interarrival  time ,  at  each 
channel  is  .32  .  This  clearly  indicates  that  the  probability  of  any 
channel  being  almost  constantly  busy  is  small  and  hence  the  probability 
of  a  call  being  lost  is  very  small. 

d)  Finally,  we  shall  consider  the  distributions  obtained  for  the 
interarrival  times  and  the  busy  times. 

With  ten  degrees  of  freedom,  the  probability  of  obtaining  a 

2  2 
X  value  greater  the  18  is  . 05  and  a  x  value  greater  the  16 

2 

is  .10.  Thus  the  x  value  obtained  for  the  interarrival  time 
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distributions,  which  are  given  in  Table  2.5  of  Appendix  II,  clearly 
indicate  the  fit  of  these  distributions  to  the  exponential  is  excellent. 
Thus,  in  this  case,  we  can  certainly  conclude  that  the  assumption  of 
calls  arriving  according  to  a  Poisson  distribution  is  reasonable. 

However,  this  is  not  the  case  when  dealing  with  the  busy  time 

2 

distributions.  With  22  degrees  of  freedom,  the  probability  of  a  X 

2 

value  greater  than  40  is  .01  and  a  x  value  greater  than  34  is 

2 

.05.  Thus,  from  the  x  values  obtained  for  the  Y  distributions, 
which  are  given  in  Table  2.5  of  Appendix  II,  we  can  only  accept  the  9 
to  10  A.M.  and  the  4  to  5  P.M.  distributions  at  the  1%  level 
of  significance.  Thus,  although  Figures  3.9  to  3.16  of  Appendix  III 
seem  to  indicate  that  these  distributions  are  tending  to  be  exponential, 
we  must  conclude  that  the  exponential  assumptions  of  the  first  two 
chapters  do  not  hold  for  our  data. 

There  are  two  main  reasons  for  this  result.  First  of  all,  we 
are  attempting  to  show  that  the  distribution  of  Y  =  T^  -  X  is 
exponential  when  we  have  no  information  concerning  the  distribution  of 
X  .  We  have  been  forced  to  assume  X  is  constant,  which  is  certainly 
not  the  case,  as  the  value  of  X  may  vary  considerably  from  call  to  call. 
However,  thirty  seconds  is  a  large  estimate  of  the  value  of  X  and  even 
in  this  case  the  fit  is  not  good. 

Secondly,  the  fact  that  these  particular  telephones  operate  on 
a  rotary  hunting  system  can  seriously  effect  the  assumptions  of 
independence.  With  this  system,  an  individual  is  aware  if  a  call  is 
waiting  for  him  on  another  line  and  thus  the  incoming  calls  can  affect 


, 

£.f  -7-. 


c 

' 

. 


53 


the  conversation  length. 

In  concluding,  we  should  keep  in  mind  that  the  behavior  and 
habits  of  the  people  who  are  using  these  telephones  cannot  always  be 
made  to  fit  the  strict  set  of  assumptions  necessary  to  carry  out  such 
an  investigation  . 


i. 
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APPENDICES 
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APPENDIX  I :  Tables  of  results  for  each 


hour  of  each  day 


1-1 


TABLE  1.1 


Observed  Frequencies  of  the  Interarrival  Times 


Between  9 

A . M .  and 

10  A. 

M.  of  Each 

Day 

Considered 

— 

Length  in 

lOths  of  a 

Second  =  r 

Number  = 

n 

r 

Oct. 

2nd 

Oct. 

3rd 

Oct . 

6  th 

Oct . 

7  th 

Oct . 

8  th 

Oct. 

9  th 

Oct. 

10th 

Oct . 

14th 

0- 

14 

17 

27 

21 

23 

16 

28 

22 

1- 

14 

19 

19 

20 

12 

10 

12 
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TABLE  1.2 

Observed  Frequencies  of  the  Interarrival  Times 
Between  10  A.M.  and  11  A.M.  of  Each  Day  Considered 
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lOths  of  a 

Second  =  r 
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TABLE  1.3 


Observed  Frequencies  of  the  Interarrival  Times 


Between  11 

A.M. 

and  12 

A.M.  of 

Each  Day 

Considered 

Length  in 

lOths  of  a 

Second  =  r 

f 

Number 

=  nr 
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2nd 
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3rd 
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Oct. 
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Oct . 
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TABLE  1.4 


Observed  Frequencies  of  the  Interarrival  Times 


Between  12 

P.M. 

and  1  P 

.M.  of 

Each  Day 

Considered 

Length  in 
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=  n 
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TABLE  1.5 

Observed  Frequencies  of  the  Interarrival  Times 
Between  1  P.M.  and  2  P.M.  of  Each  Day  Considered 
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TABLE  1.6 

Observed  Frequencies  of  the  Interarrival  Times 
Between  2  P.M.  and  3  P.M.  of  Each  Day  Considered 


Length  in 
lOths  of  a 

Second  =  r 

Number 

=  n 

r 
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Oct. 

6  th 
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Oct. 
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TABLE  1.7 


Observed  Frequencies  of  the  Interarrival  Times 


Between  3 

P . M .  and 

4  P 

.M.  of  Each 

Day 

Considered 

Length  in 
lOths  of  a 
Second  =  r 

Number  = 

n 

r 

Oct. 

2nd 

Oct. 

3rd 

Oct. 

6  th 

Oct. 

7  th 

Oct. 

8th 

Oct. 

9  th 

Oct. 

10th 

Oct. 

14th 

0- 

15 

27 

22 

24 

22 

17 

37 

1- 

16 

15 

21 

17 

18 

16 

18 

2- 

11 

16 

13 

12 

17 
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3- 
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11 
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12 

10 
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11 
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11 

7 

10 

8 

7 
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7 

7 

8 
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TABLE  1.8 

Observed  Frequencies  of  the  Interarrival  Times 
Between  4  P.M.  and  5  P.M.  of  Each  Day  Considered 


Length  in 
lOths  of  a 

Second  =  r 

Number 

=  n 

r 

1 

i 

Oct. 

2nd 

Oct. 

3rd . 

Oct. 

6th 

Oct . 

7  th 

Oct . 

8  th 

Oct . 

9  th 

Oct. 

10th 

OCt.'; 

14th 

0- 

21 

21 

16 

23 

15 

14 

26 

i 

1- 

16 

13 

15 

24 

7 

13 

21 

2- 

8 

11 

9 

9 

15 

7 

12 

3- 

10 

8 

4 

9 

8 

7 

12 

4- 

I —  . 

7 

5 

4 

9 

8 

7 

4 

5- 

1 

4 

4 

5 

11 

4 

5 

6- 
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4 

3 
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4 
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TABLE  1.9 

Observed  Frequencies  of  the  Busy  Times  Between 
9  A.M.  and  10  A.M.  of  Each  Day  Considered 


Length  in 

lOths  of  a 

Second  =  r 

Number 

=  n 

r 

1 

Oct. 

2nd 

Oct. 

3rd 

Oct. 

6  th 

Oct. 

7  th 

Oct . 

8  th 

Oct. 

9  th 

Oct . 

10th 

Oct. 

14th 

0- 

5 

2 

4 

2 

2 

3 

2 

5 

1- 

1 

5 

6 

5 

1 

2 

3 

5  ; 

2- 

5 

8 

5 

12 

11 

7 

5 

11 

3- 
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TABLE  1.9  continued 


Length  in 
lOths  of  a 

Second  =  r 

Number 

=  n 

r 

Oct . 

2nd 

Oct . 

3rd 

Oct . 

6  th 

Oct . 

7  th 

Oct. 

8  th 

Oct. 

9  th 

Oct . 

10th 

Oct . 

14th 
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TABLE  1.10 


Observed 

Frequencies  of 

the  Busy 

Times 

Between 

10  A.M. 

and 

11  A.M. 

of  Each 

Day  Considered 

(Length  in 

:10th s  of  a 

‘Second  =  r 

Number  = 

n 

r 

Oct . 

2nd 

Oct . 

3rd 

Oct. 

6  th 

Oct . 

7  th 

Oct. 

8  th 

Oct . 

9  th 

Oct . 

10th 

Oct . 

14th 

0- 

5 

1 

12 

3 

3 

1 

2 

10 

1- 

4 

7 

6 

2 

4 

5 

2 

4 

2- 

6 

7 

8 

7 

7 

9 

8 

6 

3- 

11 

9 

5 

6 

9 

6 

8 

11 

4- 

4 

4 

8 

3 

4 

2 

11 

6 

5- 

6 

8 

8 

0 

9 

5 

7 

2 

6- 

5 

10 

1 

5 

7 

2 

3 

4 

7- 

3 

3 

3 

3 

2 

2 

7 

2 

8- 

3 

5 

1 

2 

3 

2 

3 

2 

9- 

2 

3 

5 

5 

2 

0 

3 

9 

10- 
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1 
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TABLE  1.10  continued 
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Number 
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Oct. 

10th 

Oct. 
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TABLE  1.11 

Observed  Frequencies  of  the  Busy  Times  Between 
11  A.M.  and  12  P.M  of  Each  Day  Considered 


1-11 


r  ' " 

Length  in 
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Second  =  r 

Number 
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TABLE  1.11  continued 


,;Length  in 
dOths  of  a 

Second  =  r 
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Number 
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TABLE  1.12 

Observed  Frequencies  of  the  Busy  Times  Between 
12  P.M.  and  1  P.M.  of  Each  Day  Considered 


Length  in 

lOths  of  a 

Second  =  r 

Number 

=  n 
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Oct. 

2nd 

Oct. 

3rd 

Oct. 

6  th 
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Oct . 
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TABLE  1.12  continued 


Length  in 
lOths  of  a 

Number 

=  n 

r 

Second  =  r 

Oct. 

2nd 

Oct. 

3rd 

Oct . 

6  th 

Oct . 

7  th 

Oct. 

8th 

Oct . 

9  th 

Oct. 

10th 
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14th 

30- 
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0 
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TABLE  1.13 

Observed  Frequencies  of  the  Busy  Times  Between 
1  P.M.  and  2  P.M.  of  Each  Day  Considered 


Length  in 
lOths  of  a 

Second  =  r 

Number 

=  n 

r 

1 

Oct. 

2nd 

Oct . 

3rd 

Oct. 

6  th 

Oct . 

7  th 

Oct. 

8th 

Oct . 

9  th 

Oct. 

10th 

Oct .  1 
14th  j 

0- 

6 

6 

21 

3 

7 

4 

- 

x- 

6 

11 

6 

9 

6 

6 

2- 

13 

12 

13 

11 

10 

9 

3- 

10 

8 

13 

13 

9 

10 

4- 

7 

11 

9 

8 

9 

9 

5- 

8 

9 

7 

7 

8 

7 

- j 

6- 

7 

5 

6 

7 

2 

5 

— 

7- 

2 

7 

5 

4 

3 

6 

8- 

2 

2 

6 

2 

5 

1 

i 

9- 

5 

5 

4 

2 

2 

2 

i 

10- 

2 

2 

6 

0 

6 

3 

11- 

0 

3 

3 

5 

3 

2 

I 

12- 

3 

5 

2 

2 

1 

2 

I 

13- 

2 

3 

3 

3 

7 

2 

1 

14- 

3 

4 

4 

4 

1 

0 

15- 

4 

1 

5 

4 

0 

0 

16- 

t 

1 

4 

2 

0 

0 

3 

l 

1— ’ 
-J 

1 

2 

4 

1 

1 

4 

4 

H 

00 

1 

2 

2 

2 

2 

3 

1 

'  '  "I 

19- 

1 

1 

0 

0 

2 

2 

to 

o 

i 

5 

2 

0 

4 

1 

0 

21- 

1 

3 

1 

1 

1 

2 

22- 

0 

3 

1 

0 

0 

1 

23- 

0 

1 

1 

1 

3 

0 

24- 

2 

1 

1 

2 

1 

1 

25- 

0 

2 

1 

2 

0 

1 

26- 

1 

0 

1 

1 

3 

2 

27- 

4 

1 

0 

1 

1 

0 

28- 
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0 

1 

1 

1 

1 

29- 

0 

0 

1 

1 

0 

0 
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TABLE  1.13  continued 


Length  in 
lOths  of  a 

Second  =  r 

Number 

=  n 

r 

Oct. 

2nd 

Oct. 

3rd 

Oct. 

6th 

Oct . 

7  th 

Oct. 

8  th 

Oct. 

9  th 

Oct . 

10  th 

Oct . 

14  th 

u> 

o 

i 

0 

0 

0 

0 

0 

0 

i 

31- 

1 

0 

0 

1 

2 

0 

32- 

0 

0 

1 

1 

0 

0 

33- 

1 

0 

0 

0 

0 

1 

- 1 

34- 

0 

2 

0 

0 

0 

0 
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1 

0 

0 

1 

1 

0 

36- 

2 

0 

0 

0 

0 

1 

37- 

1 

0 

0 

1 

1 

0 

38- 

1 

0 

0 

0 

0 

0 

39- 

0 

1 

0 

1 

0 

1 

40- 

0 

0 

0 

0 

0 

0 

41- 

0 

0 

0 

0 

0 

1 

42- 

1 

0 

1 

0 

0 

2 

43- 

0 

0 

0 

0 

0 

0 

44- 

0 

0 

2 

0 

1 

0 

45- 

0 

0 

0 

0 

0 

0 

46- 

0 

1 

0 

0 

0 

0 

47- 

1 

1 

0 

0 

0 

48- 

0 

0 

0 

0 

0 

0 

49- 
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0 

0 

0 

50- 

0 
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0 
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51- 
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0 

0 

0 

0 
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0 

1 

0 

0 
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53- 

0 
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0 

0 

0 
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55- 
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0 
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TABLE  1.14 

Observed  Frequencies  of  the  Busy  Times  Between 
2  P.M.  and  3  P.M.  of  Each  Day  Considered 


Length  in 
lOths  of  a 

Second  =  r 

Number 

=  n 

r 

Oct . 

2nd 

Oct . 

3rd 

Oct. 

6  th 

Oct . 

7  th 

Oct. 

8  th 

Oct . 

9  th 

Oct. 

10th 

Oct . 

14th 

0- 

10 

5 

12 

4 

1 

11 

5 

| 

1- 

12 

6 

12 

10 

9 

9 

3 

2- 

15 

12 

13 

11 

9 

8 

8 

3- 

7 

8 

9 

11 

6 

3 

8 

4- 

12 

7 

7 

9 

8 

8 

7 

5- 

6 

7 

8 

7 

5 

5 

5 

6- 

4 

8 

2 

7 

2 

6 

3 

7- 

2 

5 

4 

5 

5 

3 

4 

8- 

2 

5 

4 

3 

6 

3 

1 

9- 

3 

2 

2 

3 

5 

3 

2 

10- 

1 

2 

3 

1 

3 

2 

2 

| 

11- 

4 

2 

2 

2 

2 

2 

2 

12- 

2 

3 

3 

2 

4 

4 

1 

13- 

i 

1 

2 

3 

3 

2 

3 

14- 

i 

0 

1 

5 

1 

2 

0 

15- 

i 

2 

1 

5 

1 

1 

0 

16- 

3 

2 

0 

2 

2 

0 

0 

17- 

2 

4 

3 

2 

3 

3 

0 

18- 

0 

2 

1 

2 

1 

0 

1 

19- 

0 

1 

1 

4 

1 

0 

1 

20- 

2 

0 

1 

2 

3 

1 

1 

21- 

0 

3 

0 

3 

1 

1 

0 

22- 

1 

1 

0 

1 

0 

0 

0 

23- 

1 

0 

0 

3 

0 

1 

1 

24- 

2 

1 

0 

0 

3 

1 

0 

25- 

2 

2 

0 

1 

0 

0 

0 

26- 

2 

0 

2 
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1 

1 

0 

27- 

0 

1 

0 

4 
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1 
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0 

29- 
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TABLE  1.14  continued 


Length  in 
lOths  of  a 

Second  =  r 

Number 

=  n 

r 

Oct. 

2nd 

Oct . 

3rd 

Oct . 

6  th 

Oct . 

7  th 

Oct . 

8  th 

Oct . 

9  th 

Oct . 

10th 

Oct . 

14th 

30- 

I _ 

4 

1 

0 

0 

0 

1 

0 

i  3i- 

2 

0 

0 

0 

1 

1 

0 

!  32- 

1 

2 

1 

1 

2 

1 

0 

33- 

1 

2 

0 

3 

1 

0 

1 

34- 

0 

2 

0 

0 

0 

1 

0 

35- 

0 

0 

0 

0 

0 

2 

1 

36- 

0 

1 

1 

0 

1 

0 

0 

\ 

37- 

2 

0 

0 

1 

0 

2 

0 

38- 

1 

1 

0 

0 

0 

0 

0 

1 

39- 

0 

0 

0 

0 

1 

1 

0 

40- 

0 

0 

0 

1 

2 

0 

0 

I 

41- 

0 

1 

1 

0 

1 

0 

0 

42- 

0 

1 

1 

0 

0 

0 

0 

43- 

0 

0 

1 

0 

0 

0 

0 

44- 

0 

0 

0 

0 

1 

0 

0 

45- 

0 

0 

0 

1 

1 

0 

0 

46- 

0 

1 

0 

0 

0 

0 

0 

47- 

0 

1 

0 

0 

0 

0 

1 

48- 

0 

1 

0 

0 

0 

0 

0 

1 

1 

49- 

1 

0 

1 

0 

2 

0 

0 

50- 

0 

0 

0 

0 

0 

0 

1 

I 

51- 

0 

0 

0 

0 

1 

0 

0 

I 

| 

52- 

0 

0 

0 

0 

0 

0 

0 

53- 

0 

0 

0 

1 

0 

0 

0 

54- 

0 

0 

0 

0 

0 

0 

0 

55- 

0 

0 

0 

0 

0 

0 

0 

56- 

0 

0 

0 

0 

1 

0 

0 

57- 

i  0 

0 

0 

0 
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0 

1 

58- 
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0 
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0 

0 

0 

>  59 

1  6 

2 

0 

3 

3 

4 

0 

1 

118 

110 

100 

125 

106 

95 

63 
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TABLE  1.15 


Observed  Frequencies  of  the  Busy  Times  Between 


3  P.M. 

and 

4  P.M. 

of  Each 

Day  Considered 

Length  in 
lOths  of  a 

Second  =  r 

Number 

=  n 

r 

- } 

Oct. 

2nd 

Oct . 

3rd 

Oct. 

6  th 

Oct. 

7  th 

Oct . 

8  th 

Oct . 

9  th 

Oct. 

10th 

Oct . 

14th 

0- 

2 

10 

7 

3 

8 

16 

12 

1- 

4 

4 

4 

7 

5 

2 

5 

2- 

3 

6 

15 

12 

7 

8 

7 

3- 

i 

9 

11 

4 

4 

7 

1 

8 

4- 

8 

15 

10 

9 

10 

7 

9 

fi 

5- 

6 

4 

5 

6 

3 

2 

4 

i 

6- 

6 

4 

4 

7 

3 

2 

5 

1 

L  7- 

7 

8 

3 

6 

3 

4 

2 

8- 

3 

3 

6 

5 

3 

7 

3 

9- 

6 

4 

6 

3 

5 

2 

5 

10- 

4 

3 

5 

1 

5 

2 

7 

11- 

1 

3 

4 

2 

5 

2 

3 

12- 

3 

3 

2 

1 

4 

1 

2 

13- 

1 

1 

1 

1 

2 

0 

5 

14- 

1 

4 

2 

1 

0 

1 

2 

15- 

1 

3 

1 

2 

2 

1 
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16- 
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1 

3 

1 

1 

3 
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1 

1 

2 

1 

1 

0 

1 
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1 

0 

2 

1 

0 

0 

1 

_ 1 

19- 
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0 

2 

3 

2 

3 

_ 1 
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2 

0 

3 

22- 
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2 
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2 
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2 

2 

2 
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TABLE  1.15  continued 


r - 

Length  in 

Number 

=  n 

r 

>10ths  of  a 

Second  =  r 

Oct. 

Oct. 

Oct . 

Oct . 

Oct . 

Oct. 

Oct. 

Oct. 

2nd 

3rd 

6  th 

7  th 

8  th 

9  th 

10th 

14th 

30- 

1 

0 

0 

0 

2 

1 

0 

31- 

1 

0 

0 

1 

1 

0 

0 

32- 

0 

0 

0 

2 

1 

0 

0 

33- 

0 

0 

0 

1 

0 

0 

1 

34- 

0 

0 

1 

2 

0 

2 

1 

35- 

1 

0 

1 

0 

0 

0 

1 

36- 

0 

0 

0 

1 

0 

0 

0 

37- 

0 

1 

1 

0 

0 

0 

0 

38- 

1 

1 

1 

1 

0 

0 

0 

1 

•V 

39- 

0 

0 

1 

0 

1 

0 

1 

40- 

0 

0 

0 

0 

0 

0 

0 

41- 

0 

0 

2 

0 

0 

0 

0 

42- 

0 

0 

0 

1 

0 

0 

0 

43- 

0 

0 

0 

0 

0 

0 

0 

44- 

0 

0 

0 

0 

1 

0 

0 

45- 

0 

0 

1 

0 

0 

0 

0 

46— 

1 

1 

0 

0 

0 

0 

0 

.  I 
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TABLE  1.16 

Observed  Frequencies  of  the  Busy  Times  Between 
4  P.M.  and  5  P.M.  of  Each  Day  Considered 
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TABLE  1.16  continued 
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TABLE  1.17 


Maximum  Likelihood  Estimates  of  X  for 


Each  Hour  of  Each  Day 
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TABLE  1.18 

Maximum  Likelihood  Estimates  of  ]i  for 
Each  Hour  of  Each  Day 
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APPENDIX 


II:  Tables  of  total 


hourly  results 
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TABLE  2.1 

Total  Observed  Frequencies  of  the  Interarrival  Times 
for  Each  One  Hour  Period 


Length  in 
lOths  of  a 

Number 

=  n 

r 

Second  =  r 

9-10 

10-11 

11-12 

12-1 

1-2 

2-3 

3-4 

4-5  ; 

0- 
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145 

144 

74 

190 

174 

164 

135 

i  _ 

121 

123 

124 

50 

140 

144 
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109  1 
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94 

97 

104 

66 

91 
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90 
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85 

73 

70 

40 

76 

89 

71 

58  ! 

4- 

63 

60 

48 

27 

41 

59 

56 

44 

5- 

44 

39 

34 

30 

33 

32 

43 

34  ! 

6 

33 

33 

40 

29 

31 

24 

28 
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35 

27 

28 

22 

21 

27 

25 
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19 

16 
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15 

14 

8 

12 
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13 

14 

16 

7 

2 
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10 
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3 
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TABLE  2.2 

Total  Observed  Frequencies  of  the  Busy  Times  for  Each 

One  Hour  Period 


Length  in 

Number 

=  n 

r 

 — ::  J 

lOths  of  a 

Second  =  r 

9-10 

10-11 

11-12 

12-1 

1-2 

2-3 

3-4 

4-5  1 

0- 

28 

39 

25 

16 

47 

48 

60 

34  | 
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28 

34 

33 

16 

44 

61 

31 

29  i 

2- 

64 

58 

69 

39 

68 

76 

58 
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3- 
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65 

70 

49 

63 

52 

44 
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54 

42 

59 

42 

53 

58 

68 
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45 

38 

31 

46 

43 

30 
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24 
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20 

11- 

23 
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16 

20 

14 

12“ 

23 
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TABLE  2 . 2  continued 


Length  in 
lOths  of  a 

Second  =  r 

Number 

=  n 

r 

'  "  ! 

9-10 

10-11 

11-12 

12-1 

1-2 

2-3 

3-4 

4-5 

30- 

2 

5 

1 

2 

0 

6 

4 

° 

31- 

3 

2 

0 

1 

4 

4 

3 

6  i 

32- 

7 

3 

3 

1 

2 

8 

3 

2 

33- 

5 

4 

2 

1 

2 

8 

2 

1 

34- 

5 

7 

2 

1 

2 

3 

6 

2 

35- 

2 

6 

1 

1 

3 

3 

3 

1 

36- 

4 

4 

3 

3 

3 

3 

1 

1 

37- 

5 

4 

0 

2 

3 

5 

2 

1 

38- 

1 

3 

0 

0 

1 

2 

4 

3 ! 

39- 

2 

3 

4 

0 

3 

2 

3 

3 

| 

40- 

2 

2 

0 

1 

0 

3 

0 

2  • 
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0 

1 

2 

2 

1 

3 

2 

42- 

1 

1 

3 

0 

4 

2 

1 

0  | 

43- 

1 

3 

2 

1 

0 

1 

0 
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0 

2 

1 

3 

1 

1 

0  | 

45- 

2 

2 

2 

0 

0 

2 
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TABLE  2.3 

Maximum  Likelihood  Estimates  of  A  and  y 
Each  One  Hour  Period 
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APPENDIX  III:  Graphical  representations  of 


total  hourly  results 
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FIGURE  3.4 

Interarrival  Times  -  12  P.M.  to  1  P.M. 


Graph  of  expected  frequencies 
Observed  frequencies 
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FIGURES  3.9  to  3.16:  Graphical  representations 

of  the  busy  times  for  each 


one  hour  period 
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